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Abstract 

We present improved wave functions for the ground state, Laughlin quasihole 
and quasiparticle excitations of the fractional quantum Hall effect. Lhese de- 
pend explicitly on the effective strength of Coulomb interaction and reproduce 

Laughlin's original result in the limit of no Coulomb interaction. 
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The fractional quantum Hall effect (FQHE) has attracted much interest since it was 
discovered in the 2D system of electrons subject to a high perpendicular magnetic field 
Laughlin first proposed a liquid-type ground state wave function (GSWF) 0. His wave 
functions for the ground state, quasihole and quasiparticle excitations are independent of the 
effective strength of the inter-electron Coulomb interaction. However it has been observed 
in experiments that the energy gap is greatly reduced as the effective strength of Coulomb 
interaction increases |§. 

After Laughlin's seminal work, most of the discussions on the fractional quantum Hall 
effect using a microscopic trial wave function has been confined to the strong magnetic field 
limit, at which the effect of Landau level mixing can be ignored. The only exception is the 
recent work by Price et. al. Q on the spherical geometry. Their trial wave function includes 
a term which is similar to the pseudopotential proposed by Ceperley || to include the effect 
of Coulomb interaction on the Wigner crystalization of the fermion one-component plasma. 

We present in this Letter a ground state wave function depending explictly on the 
Coulomb interaction, which is derived in a plausible manner from a Chern-Simons gauge 
field theory in the plane geometry. The modifying term is different from the term in Ref. 
||. Once a ground state wave function is given and if it is nondegenerate, wave functions 
for Laughlin quasihole and quasiparticle excitations can be written down directly follow- 
ing the Laughlin's argument of adiabatic flux insertion. To support our reasoning, we also 
have done numerical calculations. We have compared numerically the Laughlin ground state 
wave function and our trial wave function in the case of up to five particles. This calculation 
shows the superiority of our trial wave function. This positive result provides a motivation 
for further numerical efforts. 

Let's start with a non-relativistic Chern-Simons gauge field theory. The Hamiltonian is 

H = Jd 2 x (t^\D^\ 2 + + W)) , (0-1) 

where j = 1, 2, m is the effective mass of charge carrier (electron or hole), iDj = idj+aj—qAj 
and V c (ip) is the Coulomb interaction term given by 
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V e (1>) = tr\d 2 y {i>H{x) - p)—^—^H{y) - p), (0.2) 
2e J \x — y\ 

where q is the charge of the particle, e is the dielectric constant, p is the particle density of 
uniform neutralizing background charges, and tp^tp is the particle density of charge carriers 
whose mean value is p. We use the unit such that c = h = 1. In two spatial dimensions 
a spin polarized fermion can be traded with a composite of a hard-core boson and odd 
integer number of flux quanta || . We realize hard-core bosons by bosons interacting each 
other with a repulsive 5-function type interaction, which gives rise to the term (ip^ip) 2 . The 
interaction between the attached fluxes can be described by the Chern-Simons gauge field aj 
satisfying the Gauss constraint equation (l/2a)eijdi(ij = ip^ip where e±2 = —€21 = 1. Here a 
is the odd integer n times it, which is called the statistics parameter M. If B is the strength 
of external magnetic field, the gauge potential Aj = — (B/2)ejkX k in the symmetric gauge. 
Terms in the Hamiltonian are understood to be normally ordered. We will not specify the 
operator orderings here since the effects to the ground state wave function are irrelevant for 
our purpose. We emphasize that our model is not the one conventionally used in the effective 
field theory of the fractional quantum Hall effect ||. In the effective theory the hard-core 
interaction term comes from the Coulomb interaction in the microscopic theory and thus 
there is no explicit Coulomb interaction term in the effective theory ||. However in our 
derivation the hard-core interaction term turns out to have no relevance to the ground state 



wave function as shown by Haldane ||10|| . Hence we will adopt the second quantized version 
of a microscopic theory of Ezawa et.al., in which the Coulomb interaction term appears 
expilicitly JO]. Later we will take the particle limit to obtain the many-body quantum 
mechanical wave function. 

We are interested in the ground state wave fucntion describing FQHE at the filling factor 
v = 1/n, therefore p = eB /2a. It has also been known that a kind of Meissner effect occurs 



in the fractional quantum Hall liquid system [fL^l , which, in this formulation, can be stated 



as <2j — eAi = 5di = 0. It is easy to show that above two conditions are unique constant 
solution satisfying the equations of motion coming from the Hamiltonians with as well as 



without the Coulomb interaction term flTTl . 



We now take an ansatz for the boson field operator that is known to be good for the 
strongly correlated cases fT3|1 



(0.3) 



where 6{x) is the phase operator and rj(x) is the density fluctuation operator. Solving the 
Gauss constraint equation, we get 
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(x - y) 2 

Note that the statistics parameter a determines the strength of statistical interaction 
(5di(x) 5cii(y) ~ a 2 7](x) In \x — y\f](y)) , and therefore considering fermion case means that 



the interaction is strong ||14|| . This justifies our choice of ansatz. 

We substitute the ansatz in the Hamiltonian. Since fluctuations are very small in the 
liquid we keep terms only up to quadratic order of i](x) and 8(x). This is just the static ver- 



sion of harmonic approximation used in Ref . ^2] . Then the resulting quadratic Hamiltonian 
is 



H 



I 



d 2 x 



— / d 2 yr](x) In \x - y\r](y) + —r](x)r}{y) 
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+- 7T / d 2 yr](x) 
2e J \x — y 



(0.5) 



In the quadratic Hamiltonian the first three, the fourth and the last terms come from the 
kinetic term, the hard-core interaction term and the Coulomb interaction term, respectively. 
We next Fourier transform the quadratic Hamiltonian. The transfomation rule is 



/(*) = to I d xe 



f{*) 



-ikx 



/(*)■ 



(0.6) 



We here introduce the cut-off of momentum k max (= y/8np) to preserve the number of degrees 
of freedom given by the particle number in the original many body problem. We also exclude 



k = case in the summation, since this mode can be absorbed in the condensate, the 
constant part of ip'ip. The Fourier transformed Hamiltonian is then given by 



H= E \i-k 2 e{-k)9{k) 



( 1 ,9 ol vre 2 1 2 pa 2 l\ . , . . 

+ + - + — I + Vi-k)v(k) ■ 0.7 
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We note that a collective field theory, in which Fourier transforms of the density operator 

are regarded as collective fields, gives the same quadratic Hamiltonian in the large- N (N is 



the number of particles) expansion [13 



In the non-relativistic quantum field theory the field operators if>(x) and ip^(x) satisfy 
the following commutation relations: 

Mx),il>(y)] = [tf(x) 1 tf(y)] = 0, 

Mx),tf(y)] = 5(x-y). (0.8) 

We remember that i/i(x) and if)'(x) have been traded by rj(x) and 9{x) through the ansatz. 
Since we consider here the case where 7](x) is much smaller than p we may get an approximate 
representation ip(x) ~ Q ie ( x )(\ +r^/2) v /p. From these, the following commutation relations 
of r]{x) and 9{x) are obtained 

[6(x),6(y)] = [ V (x), V (y)] = 0, 

[e(x) 7V (y)] = -i5(x-y). (0.9) 

From these commutation relations we may write down directly commutation relations of 
r)(k) and 9{k) 

[9(k),9(k')] = [r ] (k),r ] (k')] = 0, 

[V(k),9(k')] = i8 k+k , (0.10) 

Due to these commutation relations the Fourier transformed Hamiltonian can be interpreted 
as the Hamiltonian of a system of 2N uncoupled oscillators with the following mass and 
natural frequency u>k for each mode 



M k = m/(pk 2 ), 
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(0.11) 



where w c (= eB/m) is the cyclotron frequency. 

Recall that a simple harmonic oscillator with Hamiltonian h = p 2 /2M + (1 /2) Mu 2 x 2 and 
commutation relations [x,x] = = [p,p], [x,p] = i, has the GSWF \l>o oc exp(— 1/2 Mux 2 ). 
Using this harmonic oscillator analogy |12| we reach at the GSWF of our system 
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(0.12) 



where Ib{= 1/VeB) is the magnetic length and A(= (e 2 /els) /& c ) is the so-called Landau- 
level mixing parameter measuring the effective strength of the Coulomb interaction. It is 



very easy to rewrite u k as uj k = oo c yl + 2v^fv\p + Aup 2 + 4z/ 2 p 4 with p = k/k max . Since 
p < 1 the apporximation u k ~ cj c [1 + v^fv\p + 0(p 2 )] converges as long as A < {y^fv)^ 1 . 
Even for the massive hole case, A is 3 ~ 5 ||[16| and thus this approximation is good even 
for v = 1/3 case. We emphasize that p M (M > 2) terms will give only contact terms of no 
physical importance. 

We now Fourier transform back to the coordinate space. 



= exp 
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At this stage the GSWF is expressed in terms of fields. To obtain the many-body quantum 
mechanical GSWF, we take the particle limit 



N 



a=l 



From Eqs. (13) and (14) we get 



(0.14) 
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(0.15) 



a<b l*a 2 M 

where z a = xJIb- We neglected the contact terms (~ 5(z a — zj,) and its derivatives) as well 
as the terms coming from the finite size effect, which break the translational symmetry of 
the original Hamiltonian. 

The ^>c is the wave function of charged boson. After (singular) gauge transformation to 
get the wave function of the original spin polarized fermions the final form can be written 
as 

1 
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where denotes the Laughlin wave function. This is the main result of our work. 

We now calculate energies for ^> L and ^> c for iV = 5. If we do not consider the particle- 
background and background-background Coulomb interactions, the Hamiltonian is 



AT N | 

- \Z a b\ 



(0.17) 
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where = d z + z/A and ■ 
be obtained analytically 
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(0.18) 
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Energies for ty L and which are defined by / d 2 z^* L ^ L (J2 a <b ^/\z a b\)/ 1 d 2 z^* L ^ L and 
/ d 2 z^* c ^ c H c / J d 2 z^* c ^ c , are calculated for N = 5, A = 2/3 case using Monte Carlo 
method, and are given by 2.05 and 1.75, respectively. That is, the energy is lowered by 
about 15%. This positive result from the five-body calculation provides a motivation for 
more numerical efforts. There is another support for the would-be superiority of our wave 
function. The total energy per particle is given in terms of the radial distribution function 
<?(|z|). Laughlin argued that the main reason why the fractional quantum Hall liquid is 
variationally superior to the Wigner crystal at v — 1/3 is the fact that g(\z\) for the liquid 



is pushed out further from the origin than <7(|;z|) for the crystal |17 |. It is evident that the 



modifying term in our trial wave function further push out <?(|z|) from the origin. Thus we 
expect our wave function will be variationally superior to the Laughlin wave function. 

Based on our improved ground state wave function, we can directly write down the 
wave functions ~ z ° an d +z ° f°r the Laughlin quasihole and quasiparticle excitations 
following the Laughlin's argument P, |T7l , |l8H of adiabatic flux insertion 



> c u = exp 



t E W 2 - ^ E rz j— 7 ] - II - z b ) 



*c + * =exp[-^El^! 2 -^E T-4-|] n ( 2 7^ - IK*- - (°- 19 ) 
One might try to derive these wave functions using the mean field method used above from a 
vortex field theory obtained via the Hubbard- Stratonovich transformation [|ll|,|l^]. However 



this is not so promising, since the linear approximation (or the harmonic approximation) is 
not good for the intrinsically nonlinear object like a vortex and moreover the vortex size is 
comparable with the magnetic length, which is a typical length of the system. 

The appearance of the FQHE at v = 1/7 has been controversal. Goldman et. al's 
experiment [|H5[] and Jain's classification |2T|] insist on the appearance. However, Lam and 
Girvin's classic work [F9| and the most updated numerical calculation show that the 



critical filling factor v c for the liquid-solid transition is larger than 1/7. In Ref. @], they 
used a trial wave function including the Landau level mixing effect, which however has a 
form different from ours in (16). It is interesting to see if our improved wavefunction gives 
the ground state energy lower than that of Ref. [f|] making v c smaller. In case v c is less than 
1/7 with our modified wave function, then the appearance of the FQHE at v = 1/7 could 
be explained theoretically. 

As discussed in the beginning of this paper, the Laughlin's wave functions for the ground 
state, quasihole and quasiparticle excitations are independent of A and therefore cannot 
explain the experimentally observed reduction of energy gap between the ground state and 
the Laughlin quasiexciton as A increases. However, our wave functions in (16) and (19) are 
explicitly dependent on A and therefore there is a chance to explain the reduction of the 
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energy gap. 

In conclusion we propose systematically derived improved trial wave functions for the 
ground state, and Laughlin quasihole and quasiparticle excitations. These include the Lan- 
dau level mixing effect explicitly. Since the validity of our wavefunction will shed new light 
on the appearance of the fractional quantum Hall effect at v = 1/7 and the reduction of 
energy gap as the effective strength of Coulomb interaction increases HJ2(1, it will be very 
interesting to continue the numerical study. 

This work was supported in part by the KOSEF and by the Korean Ministry of Educa- 
tion. 
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